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4. Statement of the Problem Studied 
Research Objectives 

We hypothesize that cavitation is a damage mechanism for blast induced traumatic brain injury (bTBI). A main objective of this project 
is to determine the conditions conducive for cavitation in cerebrospinal fluid (CSF) and corresponding tissue injury in 2-D brain slices. 
Initially, tissue surrogates will be used to develop the experimental technique, identify conditions for cavitation development, its collapse 
and associated pressures. With this knowledge, rat brain tissue will be exposed to similar conditions to develop quantitative analysis 
methods to relate neuronal damage to cavitation, pressure and strain history. Neural damage is assessed through new histological injury 
assays and mechanical property testing by indentation. With the data from applied peak pressure, rate of loading, cavitation threshold 
and histology, the ultimate goal is to develop “pressure -deformation- neural injury” maps for brain tissue. 
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Approach 


MRI scans lack resolution to identify local neural damage due to shock loading. We have developed a testing device that can resolve 
real-time incidenee of eavitation and tissue deformation in response to high rate loading in living tissues. An experimental setup to 
replicate a 2-D model of the brain has been developed. The model consists of a transparent test cell filled with artificial CSF and a brain 
(or a surrogate) sliees whieh is subjeeted to high pressure rapid loading with a polymer split Hopkinson pressure bar (PSHPB). Imaging 
data acquired with a high- speed camera is utilized to identify cavitation initiation and collapse at sub -millisecond resolution. In 
addition, the recorded tissue deformation is used to ereate strain maps for the entire region of the brain slice. In our study, we have 
concentrated our efforts on damage being caused by cavitation in the hippocampus region of the bran, which is responsible for mood 
and memory. These studies were eorrelated to damage in brain tissue slices quantified using immune-histochemical assays and seanning 
electron microscopy. 

In heterogeneous tissues such as the brain, there are limitations to using homogenous assumptions when explaining complex tissue 
interactions between different anatomical regions. There is a need for additional measures of spatial as well as transient variation in 
meehanical properties within different anatomical regions of the brain. The effect of tissue damage on ehanges in tissue stiffness is 
assessed locally in various brain regions with micro-indentation. 


Signiflcance 

Blast-induced traumatic brain injury (b-TBI) is a signature injury of wars in Iraq and Afghanistan. The root-cause of b-TBI, i.e. what 
is “broken”, in the brain during exposure to shoek loading is currently unknown. While blast waves are well known to have negative 
pressure regions, the effect of such negative pressures on brain injury has not been previously explored. In this study, we focus on 
whether incident shoek waves cause ‘cavitation’ in fluid regions surrounding the brain and subsequent neural damage in adjacent brain 
tissue. New injury mechanisms may result in novel mitigation strategies, e.g. protective gear, or identification of improved drug 
therapies specific to this injury. 


5. Summary of the Most Important Results 

We have developed new testing tools to study intracranial cavitation and the associated tissue damage eaused by cavitation. 

a. Simulated intracranial blast wave test system: A PSHPB test system was developed to generate a pressure wave that mimics a 
wave generated inside the head during a blast event, i.e. lED explosion. We were able to develop a system which can generate 
overpressure and underpressure components as well as isolate overpressure alone. 
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b. Novel controlled cavitation test system: Single bubble cavitation as may occur in the head during blast was reproduced using the 
PSHPB test system. Single bubble growth through collapse spanning over a duration of approximately 200- 600 microseconds was 
controlled and characterized over the short time scale of blast wave passage of around 5 milliseconds. Maximum bubble radii of 0.5- 
1.5 mm were generated. The ratios of overpressure to underpressure of approximately 4 to 10 times were achieved. Finally, estimated 
cavitation collapse pressures upto 124 MPa were achieved which correspond to around 125 times the measured overpressure of the 
blast wave. 

c. High resolution soft tissue response due to blast and cavitation events: It was unequivocally observed that cavitation indeed 
occurs and the collapse of this bubble indeed causes significant strain in the tissue well after the bubble has collapsed. Clear evidence 
of significant tearing and damage has been observed in brain tissue. Deformation of tissue surrogates and brain tissue slices with 
passage of the blast wave and with single bubble cavitation was measured. High speed video allowed for spatial mapping of 
deformation and strain during blast. After the bubble collapse, maximum von-Mises strains up to 0.35 were measured in a local 
shock area of 0.18 mm^. 

d. Early injury response due to blast and cavitation events: Severe neural degeneration (50 %) was detected when compared to 
control tissue after exposure to high magnitude simulated blast wave. Exposure to overpressure blast wave induced two distinct 
types of astrocytic injury responses: (a) Astrogliosis (change of quiescent astrocytes to reactive state) and (b) Clasmatodendrosis 
(autophagic degeneration of astrocytic distal processes). Preliminary scanning electron microscopy analysis revealed severe damage 
to axons within the hippocampal region following cavitation. 

e. Brain tissue mechanical properties: Brain tissues are super soft (G<1 kPa) and challenging to characterize. Baseline studies have 
been performed on tissue slices using a novel, optically-based micro-indentation test system. Viscoelastic behavior within various 
anatomical regions in rat brain tissue slices was measured and the range of instantaneous to equilibrium shear modulus was 3.8 to 
0.54 kPa in the cerebral cortex, 1.4 to 0.27 kPa in the hippocampus and 1.0 to 0.17 kPa in the caudate/putamen. 

5.1 Simulated intracranial blast wave test system 

5.1.1 Test Cell Design 

A transparent test cell with 4 mm test chamber thickness was designed to mimic the CSF layer in rat brain. A dynamic pressure sensor 
placed inside the test chamber was used to measure the pressure within the chamber during the test. The test cell incorporated designs 
features including rectangular geometry, see Figure 1 to support planar wave transit, and reduced fluid volume to restrict tissue 
displacement. 
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Figure 1. Schematic of the fluid filled test cell assembly with the piston rod and pressure sensor. 


5.1.2 PSHPB System 

Polymer split Hopkinson pressure bar (PSHPB) test system was used for generating simulated blast waves. This system is traditionally 
used to develop the stress-strain response of soft materials, e.g. foams, at high strain rates in the range of 500-5000/s. In this study, 
rather than determine constitutive response, the PSHPB device was used only to impart well-defined pressure waves to a fluid-filled 
specimen holder. In this way, the PSHPB device transferred momentum and high-rate loading to the test sample. 

A momentum trap PSHPB was employed to create a single peak shock wave with controlled intensity. The momentum trap consists of 
a flange-sleeve -rigid mass assembly as shown in Figures 2 and 3 introduced at the end of the incident bar. With this trap in place the 
generated compression wave passes through the test cell. The reflected tension wave that starts from the stress free end of the incident 
bar is absorbed by the momentum trap leading to pull back from the incident bar. 
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(a) 
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Figure 2. (a) Schematic depicting the momentum trap used in the Hopkinson bar and (b) portion 
of the incident bar with the momentum trap in PSHPB. The mass and the sleeve smoothly slide 
over the incident bar. 




















Figure 3. Momentum trap PSHPB system along with the test cell 


Using the fluid filled test cell integrated with the PSHPB system as shown in Figure 3, two distinct types of blast waves were generated. 
The first wave had similar profile as that of an open field blast wave with both overpressure and underpressure components as shown in 
Figure 4(a). The other type of wave generated was characterized by overpressure component alone, see Figure 4(b). These two wave 
profiles allowed for identification of cavitation damage in the presence of overpressure as would occur inside the head during a realistic 
blast event. 
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Figure 4. Typical simulated blast wave pressure profile in the fluid measured by the pressure sensor placed inside the chamber: (a) 
Simulated blast wave pressure profile consisting of both overpressure and under pressure components and (b) over pressure 
alone where no cavitation was observed in these experiments. 

5.2 Novel controlled cavitation test system 

For single bubble cavitation tests, the fluid-filled test chamber (Figures 1 and 3) was degassed to prevent the formation of multiple 
bubbles. To accomplish this, all parts of the test chamber were submerged in distilled water in a large beaker and degassed within a 
vacuum pump followed by soak in an ultrasonic bath to release and remove remaining small bubbles from all surfaces. A single seed 
bubble with a diameter of less than 50 pm was picked and placed at a predetermined location within the test cell. The temporal evolution 
of the bubble growth and collapse with exposure to simulated blast wave was characterized. Figure 5 illustrates a typical bubble growth 
and associated signal characteristics. 
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Figure 5. Dynamic behavior of a single bubble in response to a simulated blast wave, (a) Sequence of images showing evolution of a 
single bubble cavitation. The gas seed bubble had an initial diameter < 50 pm. (b) Corresponding pressure profile showing overpressure 
and underpressure phases. t=0 corresponds to the start of the PSHPB experiment. 
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The bubble remained approximately spherical throughout the test until collapse. The time from when the bubble started to grow until 
collapse ranged from 270 to 350 ps. The corresponding pressure profile measured within the chamber was similar to a blast shock wave 
with an overpressure peak phase followed immediately by a negative pressure phase and subsequent smaller oscillations. Minimum 
bubble size occurred at approximately the time of maximum overpressure, Po max. Maximum bubble size coincided with the minimum 
underpressure, F„ min. There was a time lag (50 to 63 ps) between true minimum pressure and maximum bubble size. This delay was 
attributed to i) the time required for the pressure wave to travel the distance between the bubble and sensor and ii) inertial delay of 
bubble growth. Bubble collapse was noted at point (n) and corresponded to disintegration of the original bubble into multiple smaller 
bubbles. 

In Figure 6, the ability to grow bubbles of various sizes and the associated pressure profiles are illustrated. Note that the maximum 
bubble size shown in the inset corresponds to the minimum negative pressure in the pressure profiles. 



Time (ms) 


Figure 6. Pressure measured in the test chamber to control growth and collapse of a single bubble. Different PSHPB air gun launch 
pressures resulted in different test chamber pressure profiles and peak negative pressures. Larger negative pressure resulted in larger 
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bubble size. The maximum bubble sizes for each test are: 1) 0.42 mm, 2) 0.97 mm, 3) 1.14 mm, 4) 1.29 mm, and 5) 1.47 mm as shown 
in the inset micrographs. 


Upon bubble collapse, a high-pressure pulse is released. Emission of a shock wave from a collapsing bubble has been well known since 
the famous theoretical work of Rayleigh (1917) on an imploding spherical cavity. The differential equation relating pressure and the 
radius of an isolated spherical bubble in an infinite, incompressible liquid is given by 

p — p 

RR + -R^=^ -^ 

2 p 

Where, R is the instantaneous bubble radius, which can be measured from cavitation bubble images captured by the high speed digital 
camera, Poo is the environmental pressure (far-field pressure) measured by the pressure transducer placed in the test chamber, and Pr is 
the pressure in the liquid at the boundary of the bubble. 

The measurable bubble collapse radius varied between 0.2 to 0.5 mm. It should be noted that bubble radius at collapse was difficult to 
determine as the bubble size approached the limit of image resolution. Therefore, it was not possible to identify the pressure at the exact 
time of collapse of the bubble from the images. To overcome this limitation, we have plotted the instantaneous bubble radius R during 
its growth and collapse as shown in Figure 7(a). A fifth-order polynomial was fit to this data from which one can calculate the 
corresponding R(t) and R(t) functions and then determine the boundary pressure Pc in the fluid for each bubble radius based on 
Rayleigh-Plesset relation as shown in Fig. 7(b). Because the radius at the instance of bubble collapse (i.e., at R=0) is unknown we fit an 
exponential curve (Pc = ke'^) for this data. The curve-fit equations for the growth and collapse of three typical bubble sizes are shown 
Figure 7(b). By setting P=0, we can now determine the elusive collapse pressure Pc at the time of bubble collapse i.e., Pc=k. 

Typical collapse pressure for a large size bubble (i.e., large negative or tensile pressures) can be as high as 100 MPa, see Table 1. For 
comparison, note from Fig.5(b) that a typical value of incident blast overpressure in our set up is only around 1 MPa and the blast under 
pressure is around 250 kPa. Thus the local pressure at the time of bubble collapse can be as high as 100 times the incident over pressure 
or around 400 times the blast under pressure. These large pressures can cause severe damage in the localized regions of the brain tissue. 
Table 1 summarizes series of tests performed and relevant blast pressures. 
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Figure 7. Bubble growth time history and calculated collapse pressure, (a) Typical bubble growth and collapse during testing. Radial 
profiles for three bubbles corresponding to low, medium and high underpressure groups are shown. 5th order polynomials were fit to 
experimental data, (b) Collapse pressure estimation method. Bubble pressure was calculated at each bubble radius for the low, medium 
and high underpressure tests shown in (a). Exponential equations were fit and pressure at r=0 was estimated to be the bubble collapse 
pressure. Pc. 


)f test parameters and bubble behavior during single bubl 

5le cavitation tests. 

Underpressure group 
(sample size) 

Low (n=10) 

Medium (n= 8) 

High (n=9) 

Pu min (kP^.) 

-86+11 

-166+39 

-315+21 

Air gun pressure (kPa) 

69 

130 

207 

P0 max (kPa.) 

851+49 

943+71 

1178+22 

Over/underpressure duration* (ps) 

586+23/157+14 

643+31/168+22 

723+24/184+13 

Lag time^ (ps) 

53+2 

58+7 

60+3 

Maximum radius of bubble (mm) 

0.47+0.12 

1.08+0.15 

1.45+0.11 

Collapse radius (mm) 

0.21+0.11 

0.24+0.13 

0.23+0.09 

Poo at collapse^ (kPa) 

6.51+3.32 

3.38+7.44 

0.51+4.65 
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Pc (10^ kPa) 

25.1+4.2 

68.6+5.3 

123.9+10.2 

P c/ P 0 max 

-29.5 

-72.7 

-125.2 


* Initial phases f Between P„ min and maximum bubble radius 
^ Chamber pressure measured away from the bubble 


5.3 High resolution soft tissue response due to blast and cavitation events 
5.3.1 Hydrogel deformation: High pressure cavitation testing 

Initial gel studies were conducted at high blast pressures on tissue surrogates. Figure 8 shows bubble initiation (b) and growth (c & d) 
near a gel edge in a degassed test chamber. These experiments showed clear evidence of shock wave generation with bubble collapse 
(f) that resulted in hydrogel indentation (g-i) well after the bubble has collapsed. 
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Figure 8. Bulk cavitation and tracking of a single bubble near a ballistic-gel edge. Figure (a) shows the incident blast wave direction 
traveling from right to left (blue arrows), (b) bubble initiation (circled), (c&d) growth, (f) collapse, and the effect on an adjacent gel 
slice (f-i). A local shock wave following bubble collapse resulted in highly localized hydrogel deformation (red arrows). Pomax =2800 
kPaand Pumin=-3S0kPa. 
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5.3.2 Single bubble cavitation on hydrogels 

Intent of the following types of tests is to measure the deformation fields when a controlled isolated bubble is allowed to grow and 
collapse in the vicinity of a hydrogel. By observing the deformation pattern utilizing the high speed imaging system, one can 
experimentally measure the strains during bubble growth and well after bubble collapse (shock wave generation) has occurred. 

Testing was conducted on a tissue phantom, speckled hydrogel slices. The test chamber was filled with distilled water, and a 1 mm thick 
hydrogel slice was placed on bottom of the test cell. It should be noted that the piston used to impart the pressure wave was in contact 
only with water in the test cell, not the hydrogel slices, throughout testing. For follow-up strain mapping, hydrogel slices were speckled 
with black stamp ink (AccuStamp, Cosco Inc) using an air brush. Methodology for the introduction and control of a single bubble was 
followed which included degassing with vacuum pump, followed by soak in ultrasonic bath and introduction of seed bubble at 
predetermined location. Next, the hydrogel slice was moved a fixed distance (0.5 mm) from the seed bubble (-0.05-0.1 mm diameter). 
Once the piston impacts the incident bar, the induced pressure wave caused the seed bubble to grow and collapse. Corresponding 
deformation of adjacent hydrogel was recorded with our high speed camera (60,000 frames per second), see Figure 10. The cavitation 
induced damage was assessed at two locations by controlling the placement of the bubble with respect to the slice; (i) bubble placed at 
the edge of the slice (Figure 9) (ii) bubble placed underneath the slice (Figure 10). 
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Figure 9. Speckled hydrogel slice deformed during single bubble growth and cavitation near the edge of the slice. 

5.3.2.1 Hydrogel deformation: Single bubble at an edge 

Using digital image correlation (DIC) and a series of high speed images, we calculated displacement of every speckle in the gel as well 
as the principle strains at each point, see Figure 10. Figure 11 shows the bubble radius during its growth and the associated principle 
strains for the nearest point in the hydrogel. Principal strains reached a peak value when the bubble grew to the maximum size and after 
the bubble collapse the gel retracted back to its original position and then due to inertia continued to move in the same direction causing 
another peak as seen in the Figure 11. The full field strain map calculated using DIC for the hydrogel revealed that the maximum strain 
occurred at material points which were in close contact with the bubble during its expansion and contraction phases. 
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As the bubble grew, the hydrogel edge was deformed by direct contact with the expanding bubble. Radial strain (srr), circumferential 
strain (see) and von Mises strain (sym) are shown in Figure 10. It should be noted that deformation and strain fields were not calculated 
in regions where the bubble overlapped the edge of the hydrogel (masked regions) since there was distortion of the observed speckle 
field due to the overlapping bubble. As these regions experience maximum strain, the values close to the edge could not be calculated. 
Overall, large deformations were captured and magnitude of calculated strains decreased with distance from the bubble center. 

The effect of local shock wave generation following bubble collapse was not seen by changes in Srr for this bubble configuration (no 
compressive strains after bubble collapse). This was because either the edge of the hydrogel was not sufficiently close to the center of 
the bubble to sense the shock wave, or the expanding bubble pushed the entire hydrogel slice away from the bubble rather than inducing 
local deformation. 
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Figure 10. Single bubble cavitation at the edge of an agarose hydrogel slice. (Top row) High speed images showing initial seed bubble 
at edge, maximum growth and collapse at selected times; (2“*^ row) radial strain, row) circumferential strain and (4* row) von Mises 
strain at the same time points. Regions at the edge were excluded due to masking by the impinging bubble. Test chamber was degassed 
before gas seed bubble (< 50 pm diameter) was placed to the side of the hydrogel (~1 mm gap distance). t=0 corresponds to the start of 
the PSHPB experiment. Simulated blast wave traveled from right to left with Pomax =1260 kPa and ,„/«=-189 kPa. 
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Figure 11. Inertia-driven changes in radial strain in the hydrogel during bubble growth and collapse, five typical experiments are shown. 
The nearest hydrogel point (in red) to the seed bubble that was not covered by growing bubble was selected. Gas seed bubbles (< 50 pm 
diameter) were placed to the side of the agarose hydrogel (~ 1 mm gap distance). t=0 corresponds to the start of the PSHPB experiment. 
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5.3.2.2 Hydrogel deformation: Single bubble underneath a tissue slice 

To better capture the effects of local shock wave generated with bubble collapse, a single seed bubble was placed underneath the 
hydrogel slice to provide more continuous contact with the tissue surrogate during the simulated blast. Figure 12 shows a series of high 
speed images for a typical experiment of bubble growth (b-c) and collapse (d-e) on the left side of each image, and the calculated strain 
field in the hydrogel at that moment on the right side of each image. 

During bubble growth, radial strains at the bubble edge are compressive. Circumferential strains at the bubble edge are tensile. Following 
collapse, highly localized compressive deformation resulted in maximum strain values {evm ranged from 0.21 to 0.51) over a short 
duration (<24 ps) at the location of the bubble center. These strains are also larger than those determined during maximum bubble 
expansion. Strains are also much larger than during passage of the compressive ID overpressure wave only (below detectable strain 
limit). Comparing with the previous test bubble configuration (edge placement. Figures 10 and 11), these localized shock strains are 
also slightly larger but over a shorter time period, see strain summary in Table 2. 
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Figure 12. Typical hydrogel deformation behavior for bubble cavitation underneath a speckled hydrogel slice. (Left image) high speed 
images showing bubble growth (b & c) and collapse (d & e) at select time intervals; (Right image) von Mises strain maps calculated 
using Die. The dark region corresponds to the bubble overlap which causes distortion in the image and so no strains were calculated in 
this region. Large compressive strains greater than 0.2 are measured due to localized shock waves at the bubble center after collapse (f). 
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Table 2. Summary of cavitation and hydrogel slice test parameters. 


Group (size) 

Edge placement 
(n=6) 

LFndemeath placement 
(n=7) 

Air gun pressure (kPa) 

207 

207 

P0 max (kPa.) 

1100+113 

1088+86 

Pumin (kPs.) 

-207+32 

-223+49 

Lag time (Ds) 

61+7 

53+6 

Maximum R (mm) 

2.35+0.34 

2.54+0.63 

Pc (kPa) 

116.3x10^+8.78x103 

129.1 x103+12.2x103 

P c/P 0 max 

105.7+1.23 

118.7+1.8 

Maximum* Srr 
with bubble expansion 

-0.15+0.03 

-0.23+0.06 

Maximum Srr 
inertially-driven tension 

0.12+0.02 

— 

Maximum* eee 
with bubble expansion 

0.06+0.01 

0.17+0.03 

Maximum* evm 
with bubble expansion 

0.17+0.05 

0.27+0.16 

Maximum evm 
after bubble collapse 

— 

0.35+0.16 

Local shock area ^ (mm^) 

— 

0.18+0.03 


* Maximum strain values in masked regions were not calculated; ^ Based on strain threshold of ^VM > 0.15 

Having characterized the single bubble test system, in the following we will present results of cavitation induced deformation and 
damage on live brain tissue slices. Figure 13 shows the effect of cavitation events on a rat brain tissue slice during a simulated blast with 
bubble initiation and growth and collapse. Note that the maximum bubble size occurs at t=1248 ps where the tissue tearing starts to 
occur well after the bubble has collapsed at time t = 1704 ps when localized tissue tearing was observed. These observations reveal the 
potential for collapse events to deform and injure soft brain tissues. 
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Figure 13. Effect of cavitation on a live rat brain tissue slice which was speckled and suspended in artificial cerebrospinal fluid. A 
simulated blast wave traveled from right to left with Po max =1045 kPa and min= -203 kPa. PSHPB testing showed bubble initiation 
(b), growth (b-d) and collapse (f). Cavitation-based damage with bubble collapse was observed in the hypothalamus region (box). The 
tissue tearing is observed well after bubble collapse (f-h). 


5.4 Early injury response due to blast and cavitation events 
5.4.1 Injury due to blast wave and cavitation: High pressure test 

The intent of this study was to establish characterize and analyze the progression of neural damage occurring during initial hours after 
exposure to high blast pressure magnitudes. The ex vivo model utilized acute rat brain tissue slices that were attached to a ballistic 
gelatin substrate and submerged in oxygenated artificial cerebrospinal fluid (aCSF). The system combined technologies of PSHPB, 
high-speed imaging, and histological methods to provide a controllable testing environment. The profile of pressure generated inside 
the chamber is graphically represented below (Figure 14). The pressure profile comprised of both underpressure and overpressure 
components and both formation, growth and collapse of bubbles for the duration of the test was observed. In this study, average peak 
fluid pressure within the test chamber reached a value of 1584±63.3 psi (Figure 14). Cavitation due to a trailing underpressure wave was 
also observed. 
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Figure 14. Artificial cerebrospinal fluid (aCSF) pressure profiles during polymer split Hopkinson pressure bar (PSHPB) testing. Profiles 
from 9 tests are overlaid. Average peak pressure was 1584.4±63.3 psi (±1 standard deviation), with an average full width half maximum 
value of 314 psec. Pressure exceeded the trigger pressure (40 psi) at 650 psec. Pressure waves were also reflected back along the incident 
bar, resulting in additional smaller impacts whose effect is assumed to be negligible compared to the F* impact pressure. The average 
peak pressure of the second wave was 483.1±30.6 psi with an average full width half maximum value of 242 psec. 

Strain maps were derived from digital video images captured by the high speed camera during PSHPB testing. Normal and shear strains 
within the brain slices were calculated at every pixel of the brain tissue images using digital image correlation method. Time-resolved 
images of tissue deformation were collected and large maximum Eigen strains (0.03-0.42), minimum Eigen strains (-0.33 to -0.03), 
maximum shear strains (0.09-0.45), and strain rates (8.4xl0^/sec) were estimated using digital image correlation (DIG), see Eigure 15. 
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Figure 15. Strain maps of the brain. The hippocampus in this is shown separately in G-I. El and E2 are the maximum and minimum 
Eigen strains, respectively. 


Neuronal injury was analyzed using a modified histological stain for degenerating neurons (Fluoro-Jade C, FJC). The protocol required 
to obtain the slices subjects the brain tissue to conditions of hypoxia, mechanical damage and injury. In order to obtain a timeline for 
mechanical testing and to determine baseline damage due to slicing protocol, effects of cell viability and tissue integrity on changes of 
mechanical properties of normal rat brain tissue slices were also investigated (Lee et al. 2012). Eigure 16 depicts the detection timeline 
and sensitivity of FJC for shorter time duration following mechanical damage. Neuronal degeneration was effectively detected after 4 
hours post slicing and tissue degradation was detected at 2 hours incubation time. 
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Figure 16. Viability testing of acute rat brain tissue slices. Hippocampal regions (DG: dentate gyrus, CAS and CAl) were tested over 2 
hour intervals (a- 0 hour, b- 2 hours, c- 4 hours, d- 6 hours and e-10 hours incubation time). Fluorescent images show bright green 
regions (red arrows) corresponding to degenerating neurons (FJC-positive). Extensive degeneration throughout the hippocampus was 
observed after 6 hours incubation. 
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Neuronal injury at 4 and 6 h after shock wave exposure in a brain tissue was quantified using Fluoro-Jade C. The total number of nuclei 
in the same region of hippocampus was also detected using DAPI staining. Quantitative cell counting methodology was applied to 
measure differences between blast exposed and control specimens between the three distinct region of the hippocampus (DG, CAl and 
CAS). The result, expressed as percentage degenerating neurons was calculated for each tissue section and sub region to assess cell type 
and tissue spatial vulnerability, see Figure 17. Statistical analyses using General Linear Models (SAS) was conducted to examine any 
statistical difference between users and between control and blast exposed samples. Neuronal injury due to PSHPB testing was found 
to be significantly greater than injury associated with the tissue slice paradigm alone (Figure 18). 



Figure 17. 4',6-Diamidino-2-phenylindole (DAPI) and Fluoro-Jade C (FJC) image analysis within CAl. (al) FJC-positive areas were 
used to identify degenerating neurons within CAl (a2) DAPI-positive areas were used to identify the total number of nuclei within 
CAl (a3) Superimposition of FJC images with DAPI images to yield the % neuronal degradation within CAl. 
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Figure 18. Neuronal degeneration in the hippocampus within control and polymer split Hopkinson pressure bar (PSHPB) test slices 
measured with Fluoro-Jade C and 4',6-diamidino-2-phenylindole (DAPI) staining. Dentate gyrus (DG), CAS, and CAl regions of the 
hippocampus were tested, and control and PSHPB were compared at 4- and 6-h incubation times. Significant differences in neural 
degeneration were measured between the control and PSHPB groups by analysis of variance (ANOVA) testing for 4 and 6h of 
incubation after testing. For the 4-h incubation time, 4 brain tissue slices were used in PSHPB tests and 2 slices were used in control 
tests, and 72 FJC images were analyzed to estimate the fraction of neuronal degeneration. For the 6-h incubation time, 8 brain tissue 
slices were tested with the PSHPB, 4 slices were used for controls, and 136 FJC images were analyzed. FJC intensities higher than 30% 
of background intensity were considered as FJC-positive neurons. Bars correspond to 1 SD (285x226 mm, 96x96 dpi). 


5.4.2 Injury due to overpressure: Low pressure test 

The overpressure component of the primary blast wave is considered to be potentially damaging to the brain. Damage to astrocytes, the 
non- neuronal cells that are implicated in neuronal metabolic maintenance, blood-brain barrier, regulation of homeostatic environment 
and tissue remodeling, via direct transient physical forces has the potential to disrupt local and global functioning of neuronal tissue. In 
this study, the hypothesis that simulated blast shock wave overpressures, transmitted through intact but isolated brain architecture, is 
capable of acute astrocyte injury at early time points was studied. The temporal and energetic profile of blast waves used in this study 
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simulated actual blast events but did not contain underpressure components. A controlled lower peak pressure (~60psi, l-2ms) was 
introduced without the characteristic negative component of the wave (Figure 19). 



Figure 19. Test chamber pressure profiles during exposure to simulated blast via PSHPB testing. Average peak pressure inside the 
chamber was 60.25 + 4.99 psi. The average full width at half maximum was 362 ps. The system allowed for a single peak loading along 
with absence of underpressure (negative pressure component). 


Based on observation from clinical symptoms implicating the hippocampus as a blast-susceptible region, injury assessment for astrocytes 
was carried out in the CAl sub-region of the hippocampus using astrocyte specific cytoskeletal protein, GFAP immunolabeling. Figure 
20, highlights the morphology of astrocytes following slicing protocol. The GFAP stain revealed normal structural features like extensive 
detailed branching along with typical star shaped cell body. 
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Figure 20. Baseline morphological features of astrocytes with GFAP labeling.Control slices were subjected to similar testing conditions 
as that of blast exposed slices without actual exposure to the blast wave or incubation. The imaging was restricted to the CAl subregion 
of both hippocampi (R= right & L = left) for a given slice. At 0 hr incubation, morphology of the astrocytes appears normal. 

Evolving astrocyte pathology emerges during the initial hours after exposure to overpressure alone, at magnitudes and durations 
comparable to real blast waves (Figure 21). Two distinct astrocyte injury profiles were identified: (a) a process similar in one respect to 
conventional astrogliosis, with an increase in GFAP immunoreactivity in GFAP-positive cells, and (b) a process comparable to 
clasmatodendrosis, an acute autolytic phenomenon, characterized by voids in cytoplasm in conjunction with rounded somata and 
disintegrating processes. The latter has not been previously associated with blast-induced injury, and neither has been previously 
demonstrated as acute responses reflecting distinct astrocyte vulnerability. 


25 





Figure 21. GFAP immunolabeled hippocampal astrocytes from control (left) and blast- exposed slices (right) as a function of post¬ 
exposure incubation time. At 1 hr. incubation time the morphological features from control (A) and blast exposed (B) look similar with 
the stellate shape and radial extension. At 2 hr. incubation qualitative visual differences can be seen between control (C) and blast 
exposed (D) astrocytes. Blast-exposed astrocytes appear hypertrophied and show enhanced GFAP staining of the cell bodies along with 
increased background staining of fine processes. After 3 hr. incubation, blast-exposed astrocytes (F) showed enhanced GFAP staining 
when compared to control slices (E). At 4 hr. incubation, clasmatodendritic features were observed in astrocytes from the blast-exposed 
slices (H) when compared to control slices (G). Astrocytes from control slices (G) showed reduced complexity of processes. Astrocytes 
from the blast-exposed slices (H) showed beaded, presumably disintegrating processes (star). Insert figures highlight are individual 
confocal planes highlighting GFAP negative voids in the cytoplasm along with rounded appearance of somata. 

Immunofluorescent GFAP distribution and intensity were analyzed across treatment groups during initial hours after the blast exposure. 
The area fraction occupied by GFAP positive cells and the intensity of the GFAP within these areas were quantified using Fiji (NIH 
Image). Contributions to overall variance were tested using a general linear model (SAS Proc GLM) for the potential effects of shock, 
incubation time, brain hemisphere, image fields (within-slice), slices (within-animal), and interactions. 
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Figure 22. Area fraction of GFAP immunoreactivity and average GFAP intensity between treatment conditions over the incubation time 
points tested. Mean + s.e.m. area fraction of GFAP immunoreactivity and average GFAP intensity as a function of incubation time. N 
within the bars indicates the slices used to collect the data. Each slice contributed to an average of 25 image fields. Main effects and 
interactions were tested in a general linear model (SAS Proc GLM). With the exception of 1 hr. incubation time, the GFAP 
immunoreactivity in astrocytes was significantly higher than control at all time points. The area fraction was found to be not affected by 
blast, time, brain side, or slice within animal. Asterisk, p <0.05 compared to control slices from same time point. 


5.4.3 Injury due to cavitation: Controlled single bubble test 

The intent of this study is to evaluate the morphological and histological changes that are induced due to cavitation events. The 
methodology developed for the controlled single bubble test is extended to study damage induced by cavitation on acute brain tissue 
slices. The pressure profile used in these test consisted of both overpressure and underpressure components (Figure 23). 
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Figure 23. Pressure profile of the input simulated blast wave with both overpressure and underpressure components. 

The region of interest was focused on the hippocampus and in these preliminary tests; the bubble was placed under the hippocampal 
region. Upon bubble growth and collapse, large deformation of one side of hippocampus was visible from the high speed images long 
after the bubble collapse (Figure 24). 
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Figure 24. Sequential images of the bubble growth and collapse, (a) start of bubble growth (b)the rectangular space indicates the bubble 
at its maximum size (c) initial bubble collapse, and (d) brain tissue deformation after large time following bubble collapse. 

Upon subsequent SEM imaging, qualitative differences were observed between the two hippocampi on the same slice. The hippocampal 
side directly influenced by the bubble (right), showed a visible surface gap in the DG layer when compared to the left side. Observing 
the surface under high magnification, long continuous processes were observed on the left side (Figure 25). At the equivalent location 
on the right side, the processes appeared to be broken and a large gap was observed. 
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Figure 25. Scanning electron microscope images of hippocampus subjected to cavitation, (a) The brain tissue slice deformation during 
the maximum bubble size, with the black rectangle highlighting the bubble placement, (b) Large field view of the right hippocampus 
that was directly impacted by the bubble growth and collapse. A fractured surface is evident in the DG layer, (c) Large field view of the 
left hippocampus. The surface appears clean without any visible damage, (d) Magnified view of the white rectangular area from image 
(a). A large gap is evident in the center, (e) Magnified view of the white rectangular area from image (b). Long continuous processes 
are visible. 


30 

















5.5 Brain tissue mechanical properties 


During blast induced TBI, brain tissue potentially undergoes complex motions such as compression, tension, shear and torsion which 
can cause axonal disconnection and distortion, necrotic cell damage and intracerebral hemorrhage. These damages ean lead to changes 
in physical structure and in turn affect the mechanical properties of the brain. However, in heterogeneous tissues like the brain, there are 
limitations to using homogenous assumptions when explaining eomplex tissue interaetions between different anatomical regions. There 
is a need for additional measures of spatial as well as transient variation in mechanical properties within different anatomical regions of 
the brain. Viseoelastic properties are useful for deseribing transient mechanieal behavior and corresponding modulus values ean vary 
with time scale. Thus, knowledge of mechanical properties of brain tissue is important in characterizing structural changes with injury 
in addition to analyzing internal stresses with deformation. 

Optical coherence tomography (OCT) is a bioimaging technique that utilizes low coherence interferometry to measure the back scattered 
light refleeted from the sample to generate cross sectional images. It is noninvasive and provides real time imaging of biological samples 
at micrometer scale. In first-of-a-kind experiments we have combined OCT with micro-indentation to measure the mechanical 
properties of brain tissue sliees. This set up allows us to test over long ereep testing times when eompared to other test systems. We 
have used this system to measure mechanical as well as optical properties within different brain regions including cerebral cortex, 
putamen, hippocampus, thalamus, and corpus callosum were measured. Indentation ereep tests determined viscoelastie properties of ex 
vivo rat brain tissue sliees in different anatomieal regions (Figure 26). 




Figure 26. Normalized creep indentation in the (a) rat cerebral cortex, (b) hippocampus, and (c) caudate/putamen over 10 min. 
Deformation depth ratio was the ratio between undeformed and deformed tissue thickness. n=6 for cerebral eortex and putamen and 
n=l for hippocampus. Data clustering was the result of diseretization introduced by limited vertieal resolution (1 pixel = 1.37pm). 
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The viscoelastic parameters and relaxation of shear modulus were estimated by fitting to an isotropic neo-Hookean model. The average 
equilibrium shear modulus for normal brain was estimated as 0.2 kPa in hippocampus region, 0.45 kPa in the cerebral cortex and 0.11 
kPa in the caudate-putamen see Figure 27. Differences in the shear modulus decay across different anatomical region in the brain can 
be attributed to the morphological and mechanical heterogeneity. Even gray regions between different anatomical regions showed 
different viscoelastic properties. This can be accounted for due to differences in size, density and structure between the anatomical 
regions. 
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Figure 27. Estimated equilibrium shear modulus in three different anatomical regions of brain tissue slices. ANOVA tests showed 
significant differences in shear moduli between the cortex, hippocampus and caudate/putamen. Statistical analysis was performed using 
SAS. For a =0.05, *** p-value <0.0002, **p-value<0.002. Bars show + ISE. 
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Effects of cell viability and tissue integrity on changes of mechanical properties of normal rat brain tissue slices were also investigated. 
Propidium Iodide (PI), FluoroJade C (FJC), DAPI and Hoechst assays were utilized to detect and changes in cell viability and 
degradation, see Figure 16 & 28. The PI stains were able to detect early necrotic cell death following tissue slicing. PI stain labeled the 
dead/ damaged neuronal cells while Hoechst stain stained all the cell bodies in a given region. Figure 28 shows a comparison of PI and 
Hoechst stains for hippocampus region of a normal rat brain. Most necrotic cell damage was detected 20 pm below tissue surface 
indicating damage due to slicing. These results provided a suitable time line for in vitro mechanical testing. 
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Figure 28. Necrotic cell death in the hippocampus of acute brain tissue slices after slicing. PI staining (red fluorescence) labeled 
damaged and dead cells. Hoechst (blue fluorescence) labeled the total cells in the same field of view. Confocal images obtained at lOx 
magnification and 10 pm depth intervals, (a) Surface of tissue slice, (b) 10 pm, (c) 20 pm, and (d) 40 pm depths below the surface of 
the tissue slice. 
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V. Technology Transfer 

Discussions were held 2010 at the Society of Biomedical Engineering (SBE) conference with Dr. Reuben Kraft on high rate induced 
deformation characteristics in gels 

Worked directly with Dr. David Moore (Defense and Veterans and Brain Injury Center) to create brain model conditions of 
clinical relevance. 

The PI gave a seminar at ARL on 28 May 2014 titled “Cavitation Induced Damage in Brain Tissue: Relevance to TBI”. Around 
40 researchers attended the seminar. Separate discussions were held with Drs. Sikhanda Satapathy and Christopher Hoppel and 
others in the group. Manuscripts of the relevant materials were also sent to ARL for distribution among interested researchers. 
Discussions were held to identify specific areas of interest to ARL so that research can continue for mutual benefit. 

Most recently, the PI met with Dr. Rohan J Banton of ARL who also expressed interest in the papers. The PI has sent those 
papers as well. 

The PI has presented his research at numerous conferences where ARL researchers such as Dr. Tusit Weerasoorya have attended 
the sessions and have interacted with the PI. His post-doc also attended the seminar at ARL and it is his understanding that she 
is working on developing similar experimental techniques for shock loading of soft tissue. 
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VII. Awards and Honors 

‘Significant Contribution ^ war J’ American Nuclear Society (ANS) - Materials Science and Technology Division (MSTD) (2014) 
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‘Member-at-Large ’ Society of Experimental Meehanics Executive Board (2015-2017) 

‘Fellow of SEM’ Soeiety of Experimental Meehanics, for “Contributions to the field of meehanics and materials through fundamental 
experimental and theoretical investigations for high strain rate, multi-axial response of novel materials, including fracture under 
extreme environmental eonditions, to inerease understanding of material behavior.” (2014) 

‘Technology Innovator Award’ University of Florida (2014). 

‘University of Florida Research Foundation Professor’ XJrmtvsXty of Florida (2013) 

‘Teacher/Scholar of the Year’ College of Engineeimg, University of Florida (2013) 

‘Researcher of the Fear ’ Meehanical and Aerospaee Engineering Department, University of Florida (2011) 

VIII. Graduate Students Involved Directly in ARO Project 

Sung, Jin Lee, Ph.D. - He has graduated 
Yu Hong, Ph.D. - He has graduated 

Saranya Canchi - She is a current Ph.D. candidate (Expected to graduate in Spring 2015) 
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Cavitation-Induced Structural and Neural Damage in Live 
Brain Tissue Slices; Relevance to TBI 


Ghatu -Suitriias-h. Univcrsi’lY of Florida. Sfi nno. f+12iinQ. HC€}. FYlfl-FYM. flRO CorePfoar-am. jSaSK 


Objective; QjuaJ^ti^y•:^avi^& 1 ^□ 4 ^ ccudirtiQfis wfttiii the 
•Eerebrcjspvnal fluid (CSF) and devdbp ' ptessure- 
defomnalian-injKJiy imaps forlbcain slices. 

Scientific Challenges: 1 . Pathogenesis of b-TlBI 
fdcMing blast expo-su re is not uei u nders-tood. 

2. Neuroimaging tedin'ques for detectiori of cavitatcii 
induced injuiy do notexbt. 

3. Experimental platfoinnsto visualkeand resolve 
r^J-time incideiice of cavitaliofi and tiss ue 
defomiatiion foiewiing Iblast exposu re are unavailable. 
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Flg.1. Ttet «■ rbr ga^eiatTg cartHUon h CS" <b) t^3ca] 
pressure pulse ard nratrun bubble stee {c] 3tia/ri a gel sJIce 
durt^ carnation ra^ltig large sbari ale'bubble celbpse 
Cega-erarg rwrons Jgireen) and toOl cells <blLie) H 
HjppocaTpus and <e) astioc.te degereiat*ori Mmrtig sbxfc 


Major AcctHn p I i s hments : 

1 . Es-tabislimentand eJiafacter^tion ofcontiolled 
single but4:4e cavitatioii fokMing blast erqposu re. 

2. Spatial mappnng of defonraation and strain on tissue 
surfogates fotowiing cavirtation. 

3. IdentiflcaliiDin of temporal prog ressbn of ineu ral and 
astrocyte ipatlioiogy fclowing blast ecposu re. 

Personnel: 2 faculty, t iresaarcli sdentEt, and 3 
graduate students (2 graduated and 1 m progr^s}. 


Army Relevance: Identifbatiion of iijuiy mechanisnis 
may iresultn novel mitlgatbti strategic-, e.g. pfotective 
gea' or ^proved drug therapies- specific to tfiisinjuiy. 
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